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Abstract — Physically-based passive microwave precipitation retrieval algorithms require a set 

24 of relationships between satellite observed brightness temperatures (TB) and the physical state 

25 of the underlying atmosphere and surface. These relationships are typically non-linear, such 

26 that inversions are ill-posed especially over variable land surfaces. In order to better understand 

27 these relationships, this work presents a theoretical analysis using brightness temperature 

28 weighting functions to quantity the percentage of the TB resulting from 

29 absorption/emission/reflection from the surface, absorption/emission/scattering by liquid and 

30 frozen hydrometeors in the cloud, the emission from atmospheric water vapor, and other 

31 contributors. The results are presented for frequencies from 10 to 874 GHz and for several 

32 individual precipitation profiles as well as for three cloud resolving model simulations of 

33 falling snow. As expected, low frequency channels (<89 GHz) respond to liquid hydrometeors 

34 and the surface, while the higher frequency channels become increasingly sensitive to ice 

35 hydrometeors and the water vapor sounding channels react to water vapor in the atmosphere. 

36 Low emissivity surfaces (water and snow-covered land) permit energy downwelling from 

37 clouds to be reflected at the surface thereby increasing the percentage of the TB resulting from 

38 the hydrometeors. The slant path at a 53° viewing angle increases the hydrometeor 

39 contributions relative to nadir viewing channels and show sensitivity to surface polarization 

40 effects. The TB percentage information presented in this paper answers questions about the 

41 relative contributions to the brightness temperatures and provides a key piece of information 

42 required to develop and improve precipitation retrievals over land surfaces. 
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43 I. Introduction 

44 

45 Because precipitation is a critical link in the Earth’s global water and energy cycles, there 

46 is an interest in retrieving precipitation (both rain and snow) in an accurate and consistent 

47 fashion on a global basis. The use of satellite observations, such as those available from the 

48 upcoming NASA Global Precipitation Measurement (GPM) mission, provides scientists with the 

49 necessary spatial and temporal coverage. Crucial for developing physically-based precipitation 

50 retrieval algorithms is an understanding of the physical relationships between satellite 

51 observations of the Earth and the state of the atmosphere and surface within the field of view. 

52 These characteristics are often interrelated, exhibiting complex spatial structures and temporal 

53 variations. Passive radiometer channels may respond to any medium within their field of view, 

54 including surface state and cloud conditions. The sensitivity to surface, temperature profiles, 

55 water vapor profiles, liquid hydrometeors, and ice hydrometeors varies depending on the channel 

56 frequency and satellite viewing angle. The lower frequency channels (10-37 GHz) are more 

57 sensitive to the surface and absorptive or emissive warming from liquid rain, while the higher 

58 frequency channels (85-874 GHz) are more sensitive to scattering from ice hydrometeors [e.g., 

59 Evans et al., 2005]. Sounding channels near the water vapor absorption lines (23, 183, 325, and 

60 448 GHz) respond primarily to the water vapor in the atmosphere and this absorption/emission 

61 continuum by water vapor generally increases with frequency. For hydrometeors, the amount, 

62 location, composition, and size distribution of the ice (and liquid) particles in the field of view 

63 produces varying brightness temperatures for different frequencies and permits retrievals of rain 

64 rates from 0.2 to 1 10 nun hr" 1 [/ Ion et al., 2008]. 

65 Precipitation retrievals from space using channels from 10-183 GHz have become a 

66 mainstay in providing rainfall rate estimates globally as evidenced by the Tropical Rainfall 
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Measuring Mission (TRMM) [e.g., Kummerow et al, 2000], the Advanced Microwave Scanning 

68 Radiometer - Earth Observing System (AMSR-E) [ Kawanishi et al, 2003], and the Advanced 

69 Microwave Sounding Units (AMSU-A and AMSU-B) [e.g., Ferraro et al, 2005]. The Ku-band 

70 (13 GHz) radar aboard TRMM [Kummerow et al., 1998] and the W-band radar aboard CloudSat 

71 [ Stephens et al, 2008] provide detailed limited-swath information about the vertical structure of 

72 cloud systems. The CloudSat radar responds primarily to cloud particles, light rain, and snow 

73 events, similar to high frequency (>90 GHz) passive radiometer channels. 

74 The TRMM and AMSR-E radiometers with channels from ~ 10-89 GHz perform 

75 particularly well over oceans when used with physically based algorithms [Lin and Hou, 2008], 

76 while retrievals over land currently use empirical approaches [Ferraro et al, 2005]. Compared to 

77 purely empirical / statistical methods, retrievals based on physical relationships are capable of 

78 providing a better understanding of the atmospheric and surface state being estimated given 

79 accurate knowledge of the necessary physical-radiometric relationships. Future retrievals over 

80 land, such as those for the NASA Global Precipitation Measurement (GPM) mission, are moving 

81 toward physically-based approaches, since these algorithms have worked well over oceans [e.g., 

82 Kummerow et al , 2009] . 

83 Recently there has also been an increased interest in estimating both ice-phase 

84 precipitation particles (snow, graupel, etc.) above the melting layer and ice particles reaching the 

85 Earth’s surface as snow in order to fully capture the atmospheric water cycle. Ice-phase 

86 precipitation detection and retrieval algorithms have been reported and shown to be useful in 

87 studying near-surface falling snow [Skofronick- Jackson et al, 2004; Ferraro et al, 2005; Chen 

88 and Staelin, 2003; Kim et al, 2008; Noli et al, 2009]. These falling snow precipitation retrieval 

89 algorithms include those that rely on neural networks, statistics, physical relationships, and/or 
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90 some combination thereof. The millimeter-wave and sub-millimeter-wave frequencies are 

91 especially sensitive both to the scattering and absorption/emission properties of ice particles, and 

92 these channels have been exploited in the above approaches. 

93 The passive signature received by satellite radiometers from falling snow events tends to 

94 be small with respect to the other contributing factors to the observed brightness temperature. 

95 However, this signature increases as the channel frequency moves from the microwave to the 

96 sub-millimeter-wave regime. When the cloud has a liquid water path (LWP) of less than about 1 

97 kg m" , the surface can contribute more than 50% of the brightness temperature for frequencies 

98 below 37 GHz. Ideally, physically-based retrievals could correct for the surface contributions 

99 and eliminate more of the uncertainty due to it. Over oceans it is much easier to determine 

100 surface emissivity and temperatures since the surface is generally relatively uniform on the 

101 typical scales of passive microwave observations. On the other hand, the emissivity and 

102 temperature of the land surface exhibits much stronger variations than ocean, and precipitation 

103 (rain or snow) at the surface can exhibit a rapid and strong influence on these quantities as snow 

104 accumulates and rain wets soil and vegetation. To determine the sensitivity to variations in the 

105 physical relationships, the contributing factors to the brightness temperatures must be understood 

106 and explained. In this paper, we perform straightforward sensitivity analyses of the associated 

107 radiometric response for multiple physical realizations of the surface and clouds (rain, anvil, and 

108 falling snow events. 

109 One approach that has been used to assist in determining the contributions from these 

110 various environmental parameters is Jacobian Decomposition [Voronovich and Gasiewski, 2004; 

111 Kim et al., 2008], With Jacobians, a tangent linear relationship is created via the use of partial 

1 12 derivatives [e.g., Voronovich and Gasiewski, 2004]. The slopes of these partial derivatives can be 
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113 positive or negative depending on whether or not the brightness temperature increases or 

114 decreases in response to changes in the atmospheric variable of interest. Other approaches to 

115 investigating the sensitivity of brightness temperatures to atmospheric and surface contributors 

116 rely on information theory or optimal estimation theory [e.g., English, 1999; Di Michele and 

117 Bauer, 2006; Bennartz and Bauer, 2003; L’Ecuyer et al., 2006]. These techniques are useful 

118 because they can be used to select the best possible instrument channels \Di Michele and Bauer, 

119 2006], to assess the impact of uncertainties [ LEcuyer et al, 2006], and to determine the 

120 sensitivity of radiances to precipitation and atmospheric profiles [Bennartz and Bauer, 2003; 

121 English, 1999]. Our infonnation content analysis approach employs the direct use of forward 

122 radiative brightness temperature (TB) calculations using the iterative or perturbation method 

123 [Lenoble, 1985] as opposed to the more commonly used adding/doubling method. 

124 The advantage of the iterative method for computing passive microwave TBs is that it 

125 provides temperature weighting functions that can be used to show the relative contributions 

126 from the combined extinction (absorption plus scattering) for the physical cloud parameters of 

127 each level in the atmosphere and for the top-of-atmosphere and surface boundaries. The 

128 temperature weighting function serves as a factor to multiply the physical atmospheric 

129 temperature in order to compute TB values [e.g., Gasiewski, 1993]. In this research, the 

130 temperature weighting function concept is exploited to analyze the TB sensitivity to the surface, 

131 atmospheric hydrometeors, the atmospheric water vapor, and other remaining constituents. In the 

132 numerical analysis, an estimate of the TB contribution (in Kelvin) from each of these cloud and 

133 surface characteristics is determined, so that the sum of these various parts equals the computed 

134 brightness temperature. The key benefit of this type of analysis is to be able to easily separate 
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135 contributions to the brightness temperatures from the physically distinct and measureable 

136 quantities in each layer of the atmosphere. 

137 This paper detennines the relationships between the key physical variables that influence 

138 passive microwave observations. We perform these analyses for selected frequencies ranging 

139 from 10 to 874 GHz, consistent with current and planned passive microwave sensors (e.g., 

140 AMSR-E, TRMM, AMSU-B, Buehler et al., [2007]). This investigation probes the TB 

141 sensitivity over several cloud classifications, such as individual profiles of convective rain, 

142 falling snow, anvil ice, as well as cloud resolving model (CRM) simulations of three snow 

143 events. Vertical atmospheric profile structures from these various cloud classifications are used 

144 as input for a robust radiative transfer model \Skofronick- Jackson et al, 2004] where nadir and 

145 53° viewed brightness temperatures are calculated. Section II describes the different case 

146 profiles, while Section III provides the procedure for computing brightness temperatures and 

147 decomposing them into quantitative percentage contributions from each source. The results for 

148 the TB and temperature weighting vectors are discussed in Section IV and Section V describes 

149 the percentage analysis results. A summary and conclusions are presented in Section VI. 

150 ii. Selection of Case Profiles 

151 While many representations of global precipitation would provide statistical results of the 

152 surface and atmospheric contributions to brightness temperature values, such results would also 

153 be diluted by the wide variety of precipitation events across the globe in all seasons and latitudes. 

154 To focus on a few key cases allows the detailed analysis needed for understanding the specifics of 

155 how surface and cloud conditions affect the resultant brightness temperature values. This 

156 information can eventually be applied in a general sense. Due to the upcoming GPM mission’s 

157 interest in the retrievals of falling snow, we emphasize snow events herein. Thus, three cloud 

158 resolving model (CRM) simulations of snow events, and six individual profiles of rain and falling 
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159 snow over both land and ocean were selected for analysis. The CRM cases allow averaging over 

160 multiple profiles of snowing events, while the individual cases allow a closer examination of 

161 specific rain, snow, and anvil cases. The CRM cases included an MM5 simulation of the 5-6 

162 March 2001 New England blizzard \Skofronick- Jackson et al., 2004; Kim et al., 2008] and a WRF 

163 simulation of the 20-22 January 2007 winter precipitation events located near the Great Lakes 

164 region [Shi et al., 2010]. Figure la shows the ice water path (IWP) image of the March 2001 

165 CRM case that deposited up to 75 cm of snow in Vermont, USA. The WRF simulation 

166 encompassed two events: (1) a shallow cloud lake effect snow on 20 January 2007 (Fig. lb) with 

167 localized snow accumulations of ~30 cm along the snow bands, and (2) a synoptic snow event on 

168 22 January 2007 (Fig. lc) with widespread snow accumulations of 5-6 cm. More complicated 

169 cases of wet/melting snowfall, multi-level clouds, and 3-D geometry effects were excluded from 

170 this analysis. 

171 The individual cases analyzed include heavy rain over ocean, light rain over land, anvil 

172 (non-precipitating) hydrometeors over ocean, light falling snow over land, moderate falling snow 

173 over land, and heavy falling snow over land. The heavy rain, light rain, and anvil cases come 

174 from an MM5 model run of Hurricane Erin [ Wu et al, 2006]. The heavy rain case was an average 

175 of 56 MM5 profiles with surface rain rates of 25±0.2 mm hr" 1 , while the light rain case was an 

176 average of 448 MM5 profiles with surface rain rates of 0.5±0.2 mm hr" 1 . The anvil case was 

177 generated from the MM5 Hurricane Erin simulation by eliminating all profiles with rain, snow, or 

178 graupel particles between 0 and 7 km, and then by selecting only profiles with 

-5 

179 I(snow+graupcl+cloud ice) above 7km to be greater than 0.12 g m" and this led to averaging 

180 over 51 profiles. The heavy snow case was determined by selecting the profile with the maximum 

181 surface snow+graupel water content (with LWP = 0) from the profiles in the MM5 cloud 
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182 simulation for the 5-6 March 2001 New England blizzard, the moderate and light snow cases are a 

183 variation of the heavy snow case as described below. 

184 The vertical profiles of hydrometeors, temperature, and relative humidity for the anvil, 

185 heavy rain, light rain, heavy snow, and light snow cases are shown in Figure 2. The heavy 

186 snowing case transitions nearly linearly with respect to height from 0.0 g m' 3 at the ~ 1 0km top of 

187 the cloud to 0.82 g in' 3 at the surface. If we assume the Marshall-Palmer [1948] particle size 

188 distribution to convert the surface snow content into a melted rain rate (e.g., RR= 19. 9M ), the 

189 RR for this profile is 15.7 mm hr" 1 or approx 6 inches an hour of a fluffy snow rate (see Table I) 

190 assuming 1 m s' 1 vertical velocity and 10-to-l ratio of fluffy snow pack to melted depth. The 

191 snowfall rate will change if the snow is considered to be more or less dense than the commonly 

192 assumed 10-to-l conversion or if the vertical velocity changes. To reduce this heavy snowfall rate 

193 into moderate and light snow rate cases, the entire vertical profile contents were multiplied by 

194 0.21 and 0.07 to obtain surface non-melted snowfall rates of 25 mm hr" 1 and 6.25 mm hr' 1 

195 respectively (Table I). 

196 In order to perform the numerical radiative transfer, the infonnation required for each case 

197 includes the cosmic background temperature, all surface inputs (surface skin temperature and 

198 emissivity), all atmospheric inputs (e.g., vertical profiles of particle size distributions, ice particle 

199 shapes and ice-air-water densities, temperature, and water vapor amounts). The simulation cases 

200 obtain most of this infonnation as part of the CRM output. Five cloud particles are allowed in the 

201 radiative transfer calculations: non-precipitating cloud water, non-precipitating cloud ice, and 

202 precipitating rain, snow, and graupel. The non-precipitating particles are monosdisperse Rayleigh 

203 particles with fixed particle size and the number density varying to account for the total 

204 hydrometeor content per layer. For the hydrometeor particle size distributions, liquid rain particles 
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205 were assumed to have spherical shapes and to follow the Marshall-Palmer [1948] particle size 

206 distribution. For the ice particles, non-spherical particles from G. Liu’s database [2004] were 

207 used. For the snow particles a relatively sparse dendritic shape was assumed (Figure 3a). For the 

208 graupel particles the denser 6-sided bullet rosette was assumed (Figure 3b). Both the dendrite 

209 snow and six-bullet graupel were considered randomly oriented with no polarized contributions. 

210 The Sekhon-Srivastava [1970] particle size distribution (PSD), N(D) = No exp (-A D), was used to 

211 obtain the A parameter in the exponential PSD. No was chosen such that ice water content (IWC) 

212 was preserved for each layer, with integration over size limited to the fixed minimum and 

213 maximum particle sizes within Liu’s non-spherical database [Liu, 2004]. 

214 The boundary conditions for the top of the atmosphere and the surface must be prescribed 

215 prior to brightness temperature computations and analysis. For all cases, the top of the atmosphere 

216 was assumed to be 2.7 K, the cosmic background temperature. The near-surface features of 

217 precipitation rate, surface temperature, surface conditions, and relative humidity for these 

218 individual cases are presented in Table I along with the liquid water path (LWP) and ice water 

219 path (IWP). For the heavy rain and anvil cases, calculations were performed over a calm ocean 

220 surface where surface emissivity could be computed using Fresnel relations [Gasiewski, 1993]. 

221 For the light rain case, a surface emissivity of 0.98 was assumed for all frequencies over forested 

222 regions following Hewison [2001], while for the individual snowing cases the surface was 

223 assumed to be deep dry snow with an emissivity of 0.82 for 10 GHz, 0.72 for 19 GHz, 0.70 for 21 

224 GHz, 0.63 for 37 GHz, 0.64 for 89 GHz, 0.72 for 166 GHz, 0.80 for 183 GHz and higher 

225 frequencies as provided by Hewison and English [ 1 999] . 

226 The surface characteristics for the MM5 and WRF models were provided by the 

227 simulations and included land surface temperature, land type, and snow depth. The land type and 
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228 snow depth of the CRM outputs was used to prescribe the emissivity of the surface. The base 

229 emissivities (£base) for CRM land type classifications of urban, crop, grass, forest, and water were 

230 classified as bare soil, close stubble, close grass, evergreen forest, and water, respectively and 

231 taken from Hewison [2001]. Because the WRF and MM5 models provide snow depth (SD) in cm 

232 as an output parameter, it was used to adjust the emissivities by the deep dry snow emissivities 

233 (£ S now) listed above for the individual snow profile cases using: 


234 


SD 

p . — p _L 

° adj ~ r, ° snow 

30 


SD^j 

30~ y 


'base 


( 1 ) 


235 Here, a 30 cm threshold is use to indicate that all low vegetation is covered by snow such that the 

236 o snow is used. For SD greater than 30 cm, only c snow is used. Since the Hewison and English 

237 emissivities were derived for localized conditions using aircraft measurements and because 

238 generally there is much variability in the emissivity, a random emissivity component of ±0.02 was 

239 added to the prescribed emissivities for each profile in the WRF and MM5 domains. These TB 

240 calculations are perfonned at the native resolution of the simulations (1km for WRF, 4km for 

241 MM5) and for each of the 207,000 and 53,000 profiles in the WRF and MM5 domains, 

242 respectively. 


243 hi. Procedure 

244 In order to detennine the percentage contributions from the surface, atmospheric hydrometeors 

245 and water vapor components of these case study clouds, a procedure to extract the separate 

246 contributions is employed. The procedure requires the use of a forward radiative transfer (FRT) 

247 model to compute brightness temperatures (TBs). During the computations, temperature 

248 weighting vectors are generated and then used to extract out the contributions from the 

249 atmospheric and surface components. 


11 



250 A. Radiative Transfer Calculations 

251 The radiative transfer equations rely on the planar-strati lied, multiple scattering based model 

252 described in \Skofronick-, Jackson et al. , 2004]. For ease in understanding the relationships, nadir 

253 brightness temperatures are the focus and presented first, followed by results at 53°. When the 

254 TBs are computed, temperature weighting functions [Gasiewski, 1993] are obtained. The 

255 brightness temperature at each frequency is the integrated sum over all heights of the product of 

256 the weighting vector value and the atmospheric temperature plus the contributions from the 

257 ground and cosmic background temperatures: 

258 ^ 

TB = f°°T(z)W(z)dz = T 0 W 0 + T CB W CB + £ T (*) W ALL (*) A 2j 

259 i= 1 (2) 

260 

261 where the subscripts 0 (at z, = 0), and CB denote ground surface and cosmic background, 

262 respectively, T(z,), W all(z,) denote the atmospheric physical temperature in Kelvin and weighting 

263 vector value in units of km' 1 , respectively, for level i of the cloud profile that consists of I total 

264 levels, and Az, is the height increment between level i and level z-1. The brightness temperature 

265 contributions (in Kelvin) from level i are the product r(z,)W |/./fz,)Az,. The weighting vectors 

266 depend on the atmospheric cloud constituents at each layer in the atmosphere. The weighting 

267 vectors W all, Wo, and W cb are defined in Gasiewski [1993] where W all is denoted W (n) . When 

268 plotted the temperature weighting functions show that the different frequencies respond to 

269 different vertical layers and cloud components (e.g., Gasiewski [1993], Skofronick- Jackson, et al., 

270 [2004]). 

271 B. Weighting Vector Decomposition and Percentage Computation 

272 The W all weighting values include the effects of every component of the vertical profile, e.g., 

273 liquid and frozen precipitating hydrometeors, non-precipitating cloud water and cloud ice, water 
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274 vapor, nitrogen, and oxygen. From Eqn. 2, the surface (S) contribution is simply ToWo, while the 

275 cosmic background (CB) contribution is TcbWcb ■ It is possible to obtain absorption-only 

276 analytical weighting function expressions (e.g., for water vapor in cloud free atmospheres as in 

277 Jackson and Gasiewski, 1995) in addition to the temperature weighting functions in Eqnl, when 

278 scattering is present analytical weighting functions for components such as precipitating 

279 hydrometeors become intractable. Therefore the approach taken here is to post-process the 

280 temperature weighting function W all (which accounts for the vertical atmospheric and cloud 

281 features) by decomposing it into separate quantitative components for studying the relative 

282 contributions from the precipitating hydrometeors (P), the non-precipitating cloud water and 

283 cloud ice (CWCI), the water vapor (RH), and other contributions (N 2 , O 2 , etc) using the following 

284 expression: 

285 W ALL (z, )T(z, ) = M P (z, y W P (z t ) + M CWCI (z f )W CWCI (z f ) + RH( Zi )W wv (z, ) + W other (z f )T{z t ) 

286 (3) 

287 where M denotes the liquid or ice water content of the subscripted variable and RH( Zi ) is the 

288 relative humidity at height z ( . This post-processing decomposition defines Wp, W wv, W cwci and 

289 W 01 her by simply renonnalizing Wall using the various extinction coefficients: 


290 

291 

292 

293 

294 

295 

296 


k-Pext (^i) T (^i) ALL (^i) 
k ext ( Zi ) M P ( Zi ) 


(K km ' 1 m 3 g 1 ) 


w ^^)= kc r ( ;- ) ^ )w r! z,) < Kkm '' 

k ext ( Z i CWCI ( Z i ) 


kt yy ( Zj ) T ( Zj ) Wall {zj) 
k K xt (z^ RH (z^ 


(K km ' 1 %■') 


(4) 

(5) 

( 6 ) 
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298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 


other ( A) — 


[kp 2 (Zj) + k N2 (Zj)} W ALL ( Zj ) 

k ext ( z i ) (km" 1 ) 


(V) 


where 


key tfei) kpext(Zi) "tk fpy(Zi)~h \i.(p\yci(Zj) ~^~\i.Q2(Zj) 4~k/v2(Zj). 


( 8 ) 


The k p ext term is the hydrometeor extinction coefficient (the sum of the absorption and 
scattering coefficients for the rain, snow, and graupel precipitating hydrometeors for layer zi), 
k cwci is the sum of the extinction coefficients for cloud water and cloud ice, k wv is the extinction 
coefficient for water vapor, and ko? and k\ 2 are the extinction coefficients for 0 2 , and N 2 . Note 
that k ext is also used non-linearly in integrals while computing W all (which is one reason 
analytical solutions for scattering weighting functions don’t yet exist). This means that Wp, 
W cwci, W wv and W ot her are minimally contaminated by the total k„, from the overall constituents, 
however, the majority of the response is captured for each weighting function component and the 
renormalization serves to elucidate the relationships between brightness temperatures and the 
underlying cloud and surface features.. 

Another point to clarify is that a large \L Pext term can be caused by a large absorption coefficient 
and/or a large scattering coefficient. In radiative transfer, large scattering coefficients typically 
result in an upwelling brightness temperature received by a satellite radiometer that is colder than 
in a non-scattering medium. The advantage of using this weighting function analysis is that it 
captures percentage contribution information from large scattering coefficients regardless of the 
reductions in brightness temperature. Figure 4 gives representative absorption and scattering 
coefficients for fixed values of hydrometeor IWC (0.01, 0.1, 0.5, and 0.75 g m" ) for five particle 
types (rain; the Liu dendrites, the Liu 6-bullet rosettes, and for comparison fluffy snow spheres 
(10% ice, 90% air) and solid snow spheres). These contents span the range of contents found in 
the vertical layers of the selected precipitation cases. Fig. 4 shows that liquid rain has a high 
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321 absorption coefficient across all frequencies, whereas the absorption coefficients for the frozen 

322 hydrometeors are 1 to 3 orders of magnitude smaller. The scattering coefficients for the various 

323 hydrometeors are similar for the 0.1, 0.5, and 0.75 g m' cases (Figs. 4b-d). For the 0.01 g in' case 

324 (Fig. 4a), the amount of ice in the frozen layer is too small to cause appreciable scattering whereas 

325 the rain does cause an order of magnitude larger scattering coefficient (but two orders of 

326 magnitude smaller than the higher content cases). Note also that the 6-bullet rosette, our “graupel” 

327 particle, has more scattering than the dendrite particle, our “snowflake” particle. One would 

328 expect that denser graupel generates more scattering than fluffy snowflake particles, and helps to 

329 justifiy our choice of these shapes for our snow and graupel particles. 

330 Once the temperature weighting vectors have been decomposed into contributions from the 

331 surface, cosmic background, atmosphere (precipitating and non-precipitating particles), and RH 

332 they can be used to obtain the value of the TB resulting from each of these components: 


333 

TB sur f- ToW'o ; TBcb ~ TcbWcb 

(9) 

334 

i 

TBp rec = T ( Zi ) W p(zi) 

i= 1 

(10) 

335 

TBcwci = X T( z i)Vy cwa(Zi) 

(ID 

336 

i 

TB rh = (z^ W wv(zi) 

i= 1 

(12) 


337 In turn, the percentages are obtained by dividing the above numbers by the total TB defined in 

338 Eqn. 2, with the remainder percentage due to oxygen, nitrogen and other components. 

339 iv. Temperature Weighing Vector and Brightness Temperature Results 

340 
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341 The temperature weighting vectors (Eq. 2, Figure 5) are used to investigate what parts of 

342 the atmosphere and boundary features contribute to the computed brightness temperature. A 

343 larger magnitude weight in W all(zi) caused by a larger magnitude extinction coefficient indicates 

344 an increased contribution from the atmosphere at those z, height levels. For example, in Figure 5 

345 for the heavy snow rate case, the magnitude of the nadir-viewed 89 GHz channel is approximately 

346 12 K km' 1 at the 1-2 km height levels, whereas 166 GHz peaks at 2-4 km with weights of ~25K 

347 km' 1 and the 183 to 874 GHz channels have larger magnitude peaks (30-40K km' 1 ), but at higher 

348 altitudes. Note that the weighting vectors are distributed across the vertical extent of the cloud 

349 indicating that the TB values are an integrated quantity over the height. 

350 When comparing the weighting functions of the heavy and light rain and anvil cases there are 

351 some similarities (e.g., lower frequencies probe to the near Earth atmospheric levels while higher 

352 frequencies are more sensitive to the upper altitudes of a cloud). There are also differences due to 

353 the amount and distribution of the hydrometeors within the clouds for the various cases. The 

354 vertical extent of the weighting function is dependent on the underlying cloud; lower LWP causes 

355 the vertical extent to be broadened (e.g., compare 89 GHz for heavy and light rain). Tracing the 

356 183±7 GHz channel weighting functions across all of the individual cases plotted in Figure 5 

357 shows that the weighting function peaks closer to the surface when the integrated quantity of 

358 hydrometeors is smaller (e.g., anvil and snow cases). This same phenomenon occurs with the 325, 

359 448, 642, and 874 GHz channels. 

360 Using the forward RT procedure, the resultant brightness temperatures for the six cases are 

361 computed and provided in Table II. This table shows the nadir brightness temperatures when 

362 each of the six profiles has the original hydrometeor profile as seen in Figure 2 and when all 

363 hydrometeors are zeroed out. The weighting function profiles and brightness temperature values 
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364 are not provided for the WRF and MM5 cases since they are composed of more than 207,000 and 

365 53,000 profdes respectively. 

366 Looking at the liquid rain cases, Table II shows that for the heavy rain over ocean case, 

367 the 10, 19, and 21 GHz channel clear-air brightness temperatures respond strongly when liquid 

368 raindrops are introduced over ocean surfaces (100K differences). The emission from the raindrops 

369 warms the TBs with respect to the radiometric ally cooler ocean surface background. There is also 

370 a strong TB cooling (—60- 1 00K) due to scattering from ice aloft in the 89 GHz and higher 

371 channels. For the light rain over land case, there is little to no sensitivity at 10-37 GHz with 

372 respect to clear air because there are few liquid hydrometeors in the lower part of the cloud and a 

373 much lower LWP (Table I). At 89 GHz and higher, again there is sensitivity and cooling due to 

374 ice scattering with 20-30K differences. The cooling due to scattering from ice at 89, 166, and the 

375 183 GHz channels is a good indicator of light rain over land. Because the anvil case has smaller 

376 ice particles at higher altitudes, only the highest channels (>325 GHz) respond with any 

377 significant change with respect to the clear air brightness temperatures. 

378 The TB values associated with the snowing cases in Table II show that for frequencies 

379 between 10 and 89 GHz, there is less than a 15K difference between the heavy, moderate, and 

380 light snow, and clear air brightness temperature values. This means that these channels are nearly 

381 insensitive to falling snow at these rates and for these cloud profiles. As the frequency increases, 

382 there is greater sensitivity to falling snow particles: At 166 GHz, there is a 48K difference 

383 between clear air and the heavy snow case, but the differences are insignificant between clear air 

384 and the moderate and light snow cases. The 183±7 GHz channel are sensitive (~35K) to the deep 

385 (10km tall) heavy snow clouds. The 325, 642, and 874 GHz channels show strong sensitivity to 

386 falling snow (e.g., at 325 GHz the TB show a nearly 50K differential between the heavy snow and 
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387 clear air cases). Between the clear air case and the light snow case for all channels the TB 

388 differences are less than 5K and less than 10K between the moderate and clear air cases. 

389 Given the low sensitivity for the moderate snow and especially the light snow cases for 

390 current satellite channels, thresholds of retrieval detection and estimation will be dependent on the 

391 snowing cloud being observed and available ancillary data. Also, since during retrievals one does 

392 not know the true state of the surface, nor the cloud microphysics, there are additional unknowns 

393 that can cause errors. Specifically, at 89 GHz, each 0.1 increment in surface emissivity error 

394 causes a 25 K change in brightness temperature for clear air profiles (Figure 6). Similarly, for a 

395 fixed surface emissivity, an error of 10°C for specifying the land surface skin temperature 

396 (effective emitting temperature) can cause a 2-3 K change in the TB values. As expected, Fig. 6 

397 shows that the surface characterization has a large impact on brightness temperature values, 

398 especially for clear-air or nearly clear-air conditions for window channels at 10, 37, 89 and 166 

399 GHz. The figure also shows that the 183 GHz channels are less sensitive to surface features due to 

400 their sensitivity to water vapor in the profile. 

401 v. Percentage Analysis 

402 By piecing together the contributions from the precipitating hydrometeors within the 

403 cloud, the atmospheric water vapor and the surface, we can have a better understanding of how all 

404 the components work together to produce a brightness temperature related to the stonn scenario. 

405 Eqns. 10, 11, and 12 provide the basis for determining the percentages for the six individual and 3 

406 simulated cases as reported in Figures 7-10. The percentage of the brightness temperature 

407 resulting from the cosmic background for the individual cases is not shown in the figures, but is 

408 no more than 1.3% (at 10 GHz for the anvil case) and is closer to 0.0-0. 3% for the other channels. 
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409 Figure 7 shows the percentage contributions for the heavy rain, light rain, and anvil cases. 

410 For the heavy rain case (Fig. 7a) the hydrometeors (liquid rain and frozen hydrometeors) provide 

411 more than 50% of the brightness temperature value. For channels at 183±1 and 448±1 GHz, the 

412 water vapor is a strong contributor to the brightness temperatures. It is only at 10 GHz where the 

413 surface significantly contributes (-35%) to the brightness temperature. For the light rain case (Fig. 

414 7b), the hydrometeors contribute less than for the heavy rain case (because their integrated 

415 quantity is less than that of the heavy rain). As such, the surface contributes more than 70% for all 

416 frequencies less than 89 GHz. We also note that there is an increasing contribution from the non- 

417 precipitating cloud water and cloud ice as the channel frequency increases. 

418 To explore the differences between heavy rain over ocean versus over land, the FRT 

419 calculations were performed with the heavy rain case over the same land background (forested) as 

420 used for the light rain case. When these percentage contributions were compared to the heavy rain 

421 over ocean, the most distinct change was at 10 GHz, where the contribution from the 

422 hydrometeors was greater for the ocean case than the land heavy rain case (33% hydrometeor 

423 contribution from land case versus a 62% hydrometeor contribution from the ocean case). This is 

424 due to the reflective surface of water (lower emissivities) which pennits energy downwelling from 

425 hydrometeors in clouds to be reflected at the surface and included in the brightness temperature as 

426 measured from an aircraft or satellite. This is an important finding since snow covered surfaces 

427 also have more reflective surfaces, and this fact will need to be accounted for in any FRT of 

428 retrievals of precipitation over snow-covered surfaces. 

429 For the anvil case in Fig. 7c, at low frequencies that are typically sensitive to the liquid 

430 drops (which don’t exist in the anvil case) the brightness temperatures are more sensitive to the 

431 surface. At higher frequencies the anvil percentages due to the water vapor in the profile are 50- 
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432 90% and at the highest frequencies (642 and 874 GHz) the percentages due to the hydrometeors 

433 are nearly 10-40% because these channels are especially sensitive to small ice particles as seen in 

434 anvil clouds [ Evans et al., 1998; Evans et al., 2005], 

435 In Table III, the percentage contributions for the surface, hydrometeors and relative 

436 humidity are provided for radiative transfer at a 53° viewing angle. Note that these percentages do 

437 not add to 100% because the contributions from CWCI, 02, N2, and cosmic background are not 

438 presented. In this table, the largest differences at 53° are at the lower frequencies where the 

439 oceanic surface polarizations affect the TB values at vertical (V) and horizontal (H) polarizations, 

440 while at higher frequencies the surface polarization percentages do not change because these 

441 channels are not sensitive to the surface. The V and H polarization percentages of the light rain 

442 case are equivalent because there is no surface polarization for land. Furthermore the non- 

443 spherical snow and graupel are randomly oriented such that they do not contribute to polarized 

444 TB values. Comparing these percentages to those in Figure 7 (for nadir viewing) there is a slight 

445 increase in the percentages from the precipitating hydrometeors and RH for the V & H 

446 percentages because of an increase in slant path for the 53° viewing angle calculations. 

447 For Figure 8, which shows the percentage contributions for the heavy, moderate, and light 

448 snow cases, there is a progression of decreasing percentages from the hydrometeors as the 

449 snowing rates decrease. Figure 8a, shows the percentages due to hydrometeors range from 20 to 

450 90% for the heavy snow case (channels 89 to 874 GHz). For the moderate snow case (Fig. 8b) at 

451 these frequencies, the hydrometeor percentages are between 5 and 30%, whereas for the light 

452 snow case (Fig. 8c), the percentages from the hydrometeors drop to 1 to 10%. There is a 

453 corresponding increase in the percentages from water vapor in the cloud. Not surprisingly, as 

454 shown in Fig. 8 for the snowing profiles (non-raining cases), the contributions for channels 10 to 
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455 89 GHz show a significant percentage from the surface. Figure 4 gives a key piece of information 

456 as to why there is less of a contribution from the hydrometeors (and more from the land) when 

457 only frozen hydrometeors are present in the profile. Figure 4 shows that the absorption and 

458 scattering of rain is higher than that of frozen particles for frequencies up to 100 GHz. Thus when 

459 liquid particles are present, the K Pext tenn in Eqn. 4 is high and hence the value of TB prec in Eqn. 

460 1 1 is elevated. 

461 Table IV provides the percentages for the snowing cases. The nadir-viewing angle 

462 percentages are those that appear in Fig. 8. These nadir viewed values are compared to the 53° 

463 angle percentages. Since the land surface component and randomly oriented snow particles have 

464 no polarization dependence, the V and H percentages are not different and only one is shown in 

465 Table IV. The 53° viewing angle results show a slight increase in percentages from the 

466 hydrometeors with corresponding deceases in surface percentages and RH percentages. This 

467 slight increase in percentages from the hydrometeors is due to the increased slant path at 53°. 

468 Figure 9 provides the percentages for the WRF 20 January 2007 lake effect snow case 

469 (panel a) and the 22 January 2007 synoptic snow case (panel b). In order to partition the 

470 percentage results for the 207,000 profiles within each of the WRF CRM simulation, an average 

471 was taken only over those profiles with IWP of 0.5 to 1.0 kg m' (18,160 profiles for 20 January 

472 2007 and 71,420 profiles for 22 January 2007). We can then intercompare the lake effect case to 

473 the synoptic case to illuminate any differences in percentage based on the structure of the vertical 

474 profiles. Thus comparing Fig. 9a (lake effect) to Fig. 9b (synopotic), the most noticeable 

475 differences are for channels 166, 183, and 325 GHz where the synoptic snow event case is less 

476 sensitive to the surface and more sensitive to the snow particles. This is likely because the lake 

477 effect events are shallow (3 to 4 km) clouds, whereas the synoptic events are 6 to 7 km tall. Thus, 
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478 despite both plots being averaged over the 0.5 to 1.0 kg m' IWPs, the taller clouds impact the 

479 weighting functions over a larger vertical extent than the shallower clouds. 

480 Finally, in Figure 10 the percentages from the MM5 heavy blizzard case are shown. For 

481 this case, percentages for an IWP of 0.5 to 1.0 kg m' (5325 profiles averaged, Fig. 10a) and for 

482 an IWP of 9.5 to 10.0 kg m' 2 (5 profiles averaged, Fig. 10b) were plotted. At low frequencies, the 

483 dominating factor is evident by the high percentages from the surface for all cases shown in 

484 Figure 10 (i.e., there is no liquid water to affect low frequency TB values). Since this is a deep 

485 snowing cloud (to ~10 km altitude) there is both a greater effect from water vapor in the cloud 

486 (for those water vapor channels of 183, 325, and 448 GHz) and a greater effect from 

487 hydrometeors in the profile. Thus there are competing factors that modify the percentages from 

488 the snow in the profile and the water vapor in the profile. A detailed analysis to determine the 

489 thresholds of detection for falling snow cases needs to be undertaken to link TB depressions to the 

490 underlying falling snow event characteristics (e.g., cloud top, IWP, particle sizes and shapes, 

491 water vapor profile, etc.). 

492 vi. Conclusions 

493 This work has provided an analysis method to determine the contributions to brightness 

494 temperature values from surface emission, from atmospheric cloud hydrometeors, from cloud 

495 water plus cloud ice, from the water vapor within the vertical profile, from the cosmic 

496 background, and from CF and No. This analysis was performed for nadir-viewing and 53° viewed 

497 observation channels spanning 10 to 874 GHz, including both window and water vapor sounding 

498 channels. These channels are relevant to current and confirmed future missions as well as 

499 channels yet to be approved for a satellite mission. Figs 7-10 show the percentage of the 

500 brightness temperature values resulting from the various contributing factors for several 
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individual profile cases (heavy rain over ocean, light rain over land, anvil over ocean, heavy snow 
over land, moderate snow over land, and light snow over land) as well as for the average 
percentages resulting from three cloud resolving model simulations (for a lake effect snow event, 
a synoptic snow event, and the 2001 New England blizzard). Across all these cases several key 
points became evident: 

1) The low frequency channels (<89 GHz) are extremely sensitive to the surface (with > 40% 
contribution to the TBs) except in the presence of high LWP (e.g., heavy rain cases). 

2) The sounding channels have a large contribution from the water vapor in the vertical 
profile, when there is a limited quantity of hydrometeors in the profile. 

3) The higher frequency channels (both sounder and window channels) are sensitive to the 
frozen hydrometeors in the cloud. 

4) Low emissivity surfaces (water and snow-covered land) pennit energy downwelling from 
hydrometeors in clouds to be reflected at the surface and therefore increase the percentage 
of the brightness temperature resulting from the hydrometeors. 

5) The type and vertical structure of a snow event (e.g., lake effect, synoptic, blizzard) affects 
the percentages resulting from the surface and hydrometeors. 

6) The slant path at a 53° viewing angle increases the percentage contributions from 
hydrometeors relative to nadir viewing channels and the percentages from the surface show 
sensitivity to surface polarization effects. 

This work provides quantitative percentage values of the sources of brightness temperature 
signals and corroborates earlier investigations indicating that low frequency channels (<90 GHz) 
are sensitive to the surface and liquid particles, and that the high frequency channels are sensitive 
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523 to frozen particles and to atmospheric water vapor. These relationships are complicated, non- 

524 linear, and vary for each channel and remotely sensed precipitation event. The number of degrees 

525 of freedom in the physical state of precipitation events is larger than the degrees of freedom 

526 available from remotely sensed observations (for current and near future platforms). Thus for 

527 retrievals over land that are physically based, ancillary data will be extremely useful in reducing 

528 the unknowns and assumptions. Primarily, infonnation about the surface emission as well as the 

529 vertical temperature and water vapor profiles are needed. It is probable that numerical weather 

530 prediction models can provide adequate vertical information. On the other hand, surface emission 

531 may require concentrated research and modeling efforts especially over certain surfaces (e.g., 

532 snow covered surfaces and deserts). While awaiting globally accurate models of surface emission, 

533 other approaches might include statistical analysis of climatological emissivities, retrievals of 

534 emissivity during clear air overpasses (with dynamic modeling of the surface as it rains or snows), 

535 or even assuming that clear air brightness temperatures just prior to a precipitation event represent 

536 the contributions from these environmental conditions. Since surface emission plays a large role 

537 for the typical frequencies used for precipitation sensors (10 to 183 GHz), future work should 

538 focus on obtaining realistic and validated surface emission globally as concurred by Noh et al., 

539 [2009]. Nevertheless, the results presented herein are useful for retrieval algorithm development 

540 for upcoming missions such as GPM, especially over land surfaces and for falling snow events. 
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636 Figure Captions: 

637 

638 Figure 1 : The IWP in kg m' 2 for the (a) MM5 March 2001 blizzard, (b) the 20 January 2007 lake 

639 effect snow stonn, and (c) the 22 January 2007 synoptic snow event. 

640 

641 Figure 2: Vertical profdes of hydrometeor liquid and ice water contents (upper plots) and the 

642 temperature and relative humidity (lower plots) for the anvil, heavy rain, light rain, heavy snow 

643 and light snow cases. 

644 

645 Figure 3: For this work, the (a) dendrite particle shape used for snow particles, and (b) the six- 

646 bullet rosette shape used for graupel particles. The non-precipitating cloud ice is assumed 

647 spherical. (From Liu, 2004) 

648 

649 Figure 4: Coefficients for absorption (solid lines) and scattering (dashed lines) for liquid rain, 

650 dentritic snow, six-bullet rosettes, solid ice spheres, and fluffy ice spheres (10% ice, 90% air) for 

651 hydrometeor contents of (a) 0.01 g m' 3 , (b) 0.1 g m' 3 , (c) 0.5 g m' 3 , and (d) 0.75 g m' 3 . 

652 

653 Figure 5: Temperature weighting vectors Wall( z <)T( z <) for the nadir-viewed anvil, heavy rain, 

654 light rain, heavy snow, and light snow cases. Upper panels are for TRMM frequencies (10-89 

655 GFIz), middle panels are for AMSU-B frequencies (89-183±7 GHz), and lower panels for the 

656 submillimeter-wave channels (325-874 GHz). 

657 Figure 6: Brightness temperature values using a standard clear-air winter atmosphere profile for 

658 varying surface emissivities and surface temperatures at (a) 10 GHz, (b) 37 GHz, (c) 89 GHz, (d) 

659 166 GHz, (e) 183±3 GHz, and (f) 183±7 GHz. The X-axis indicates the values used for surface 

660 emissivity, while the colored lines indicate surface temperature. 

661 Figure 7: Percentage contributions from the surface (green), the hydrometeors (red), the relative 

662 humidity (light blue), the cosmic background (<1% for all channels except 10 GHz where it is 

663 1.3%)+cloud water+cloud ice (purple), and the atmospheric oxygen and nitrogen (gray) to the 

664 resultant nadir brightness temperatures for (a) the heavy rain over ocean case, (b) the light rain 

665 over land case, and (c) the anvil profde case. 

666 Figure 8: Similar to Fig. 7, percentage contributions for the single profde (a) heavy snow, (b) 

667 moderate snow, and (c) light snow cases. 

668 

2 

669 Figure 9: Similar to Fig. 7, percentage contributions averaged over IWP of 0.5 to 1.0 kg m'~ for 

670 the (a) WRF 20 January 2007 lake effect case and (b) the WRF for 22 January 2007. 

671 

672 Figure 10: Similar to Fig. 7, percentage contributions of a March 2001 snow blizzard MM5 

673 simulation averaged over (a) IWP of 0.5 to 1 .0 kg in' 2 and (b) IWP of 9.5 to 10.0 kg m' 2 . 

674 

675 

676 

677 

678 

679 

680 
681 
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TABLE I. 


Surface Conditions for the Case Profiles. 


Profile 

Type 

Conditions 

T 

(°C) 

RH 

(%) 

Liquid 

Surface 

Rate 

(mm hr' 1 ) 

LWP 
(kg in 2 ) 

IWP 
(kg in' 2 ) 

Heavy Rain 

Ocean 

Calm Winds 

22.5 

100.0 

25.0 

5.37 

1.77 

Light Rain 

Land 

Forested 

23.1 

95.2 

0.5 

0.49 

0.83 

Anvil 

Ocean 

Calm Winds 

23.1 

88.9 

0.0 

0.0 

0.28 

Light Snow 

Land 

Deep dry snow 

-6 

93.8 

0.5 

0.0 

0.32 

Moderate Snow 

Land 

Deep dry snow 

-6 

93.8 

2.5 

0.0 

0.94 

Heavy Snow 

Land 

Deep dry snow 

-6 

93.8 

15.7 

0.0 

4.24 


683 

684 

685 

686 

687 

688 table ii. Brightness Temperatures For the Heavy, Moderate, and Light Snow Cases 

689 AND THE HEAVY RAIN, LIGHT RAIN, AND ANVIL PROFILES FOR “ALL” HYDROMETEORS INCLUDED 

690 AND CLEAR AIR “C A” CONDITIONS . 


Freq. 

GHz 

Heavy Snow 
Over Land 

Moderate 

Snow 

Over 

Land 

Light 

Snow 

Over 

Land 

Heavy Rain 
Over Ocean 

Light Rain 
Over Land 

Anvil 

Over Ocean 

Case 

AH 

CA 

All 

All 

All 

CA 

All 

CA 

All 

CA 

10 

217 . 1 

216.7 

216.8 

222.4 

212.0 

121 . 9 

286.5 

286.5 

121.4 

121.4 

19 

189.0 

188.5 

188.4 

189.4 

266.1 

165.0 

295.1 

296.1 

159.3 

159.6 

21 

196.1 

196.1 

197.7 

198.0 

268.2 

205.4 

289.6 

290.2 

195.5 

196.4 

37 

184 . 8 

178.6 

181 . 8 

175.2 

250.0 

178.2 

289.5 

294 . 1 

173.1 

173.1 

89 

216.2 

199.0 

207.7 

199.2 

192.7 

255.5 

271.5 

291 . 9 

247 . 6 

247.7 

166 

206.0 

252.2 

247.4 

250.9 

174.7 

278.5 

250.2 

279.6 

278.7 

280.4 

183±1 

240.8 

242 . 8 

243.0 

242.0 

232 . 6 

242.5 

240.6 

242.5 

241.3 

243.1 

183±3 

239.3 

252.3 

250.0 

253.7 

223.4 

254 . 6 

249.2 

254 . 8 

256.5 

255.5 

183±7 

225.0 

260.9 

253.5 

256.2 

196.0 

266.7 

250.6 

267.2 

267.5 

268.4 

325±7.0 

211 . 8 

258.4 

242.7 

253.0 

191.2 

262.2 

229.6 

262.7 

256.8 

263.8 

325±9.5 

208.1 

259.9 

241.2 

252.0 

189.0 

264.4 

227.0 

264 . 9 

258.5 

266.1 

448±1.5 

229.1 

231 . 6 

223.3 

230.9 

216.9 

232.3 

227.4 

232 . 1 

232 . 1 

233.1 

448±7.2 

226.92 

244.4 

236.4 

240.7 

203.4 

247.2 

232.4 

247.2 

241 . 6 

247 . 6 

642 

227 . 6 

245.0 

237.7 

235.3 

206.4 

248.2 

229.3 

248.2 

237.7 

248.6 

874 

230.0 

246.3 

237.4 

239.9 

210.3 

249.9 

227 . 6 

249.9 

234 . 1 

250.4 
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692 

693 

694 


table hi. Percentages for the Heavy Rain, Light Rain, and Anvil cases at a 53 ° 

VIEWING ANGLE. NOTE THAT THE PERCENTAGES FOR V AND H POLARIZATION FOR LIGHT RAIN 

over Land are the same since the land surface is not polarized. 


Freq. 

[GHz] 

Heavy Rain 
over Ocean 
V-pol 

Heavy Rain 
over Ocean 
H-pol 

Light Rain 
over Land 
V-pol 

Light Rain 
over Land 
H-pol 

Anvil over 
Ocean 
V-pol 

Anvil over 
Ocean 
H-pol 

Source 

surf 

hydr 

RH 

surf 

hydr 

RH 

surf 

hydr 

RH 

surf 

hydr 

RH 

surf 

hydr 

RH 

surf 

hydr 

RH 

10 

33 

64 

2 

16 

81 

2 

96 

1 

2 

96 

1 

2 

93 

0 

4 

82 

0 

10 

19 

5 

87 

8 

2 

89 

8 

76 

5 

17 

76 

5 

17 

66 

0 

31 

41 

0 

53 

21 

2 

79 

18 

1 

80 

18 

56 

6 

37 

56 

6 

37 

45 

0 

53 

23 

0 

73 

37 

0.1 

96 

2 

0.1 

96 

2 

64 

19 

13 

64 

19 

13 

67 

0 

24 

43 

0 

42 

89 

0 

95 

2 

0 

95 

2 

14 

56 

26 

14 

56 

26 

34 

0 

59 

18 

1 

73 

166 

0 

93 

4 

0 

93 

4 

0.3 

63 

35 

0.3 

63 

35 

1 

3 

94 

1 

3 

94 

183±1 

0 

51 

47 

0 

51 

47 

0 

11 

86 

0 

11 

86 

0 

3 

95 

0 

3 

95 

183±3 

0 

78 

20 

0 

78 

20 

0 

35 

62 

0 

35 

62 

0 

5 

92 

0 

5 

92 

183±7 

0 

89 

8 

0 

89 

8 

0 

55 

42 

0 

55 

42 

0 

5 

92 

0 

5 

92 

325±7.0 

0 

87 

7 

0 

87 

7 

0 

75 

19 

0 

75 

19 

0 

17 

78 

0 

17 

78 

325±9.5 

0 

87 

6 

0 

87 

6 

0 

77 

18 

0 

77 

18 

0 

16 

79 

0 

16 

79 

448±1.5 

0 

48 

47 

0 

48 

47 

0 

14 

82 

0 

14 

82 

0 

8 

88 

0 

8 

88 

448±7.2 

0 

68 

18 

0 

68 

18 

0 

53 

38 

0 

53 

38 

0 

28 

62 

0 

28 

62 

642 

0 

68 

13 

0 

68 

13 

0 

58 

28 

0 

58 

28 

0 

38 

48 

0 

38 

48 

874 

0 

69 

7 

0 

69 

7 

0 

59 

19 

0 

59 

19 

0 

42 

34 

0 

42 

34 























































































































































































































































































































696 


697 table iv. Percentages for the Heavy, Moderate, and Light Snow cases at a Nadir and 

698 53° VIEWING ANGLE. NOTE THAT THE PERCENTAGES FOR V AND H POLARIZATION ARE THE SAME. 



Nadir Viewing Angle 

53° viewing angle 

Freq. 

[GHz] 

Heavy Snow 
over Land 

Moderate 
Snow over 
Land 

Light Snow 
over Land 

Heavy Snow 
over Land 

Moderate 
Snow over 
Land 

Light Snow 
over Land 

Source 

surf 

hydr 

RH 

surf 

hydr 

RH 

surf 

hydro 

RH 

surf 

hydr 

RH 

surf 

hydr 

RH 

surf 

hydr 

RH 

10 

97 

2 

0 

98 

0 

0 

98 

0 

0 

95 

3 

0 

97 

0 

0 

97 

0 

0 

19 

91 

5 

3 

92 

0 

5 

93 

0 

5 

86 

8 

4 

88 

1 

7 

88 

1 

7 

21 

84 

7 

8 

86 

1 

11 

86 

0 

11 

76 

10 

12 

79 

1 

17 

79 

1 

17 

37 

84 

11 

1 

87 

2 

4 

86 

2 

4 

77 

16 

2 

80 

3 

6 

79 

2 

6 

89 

58 

37 

2 

71 

8 

13 

71 

6 

15 

45 

48 

3 

59 

11 

19 

59 

9 

21 

166 

13 

80 

6 

41 

13 

44 

42 

8 

48 

10 

82 

6 

25 

16 

56 

26 

10 

61 

183±1 

0 

36 

63 

0 

1 

98 

0 

0 

99 

0 

31 

68 

0 

1 

97 

0 

0 

98 

183±3 

0 

63 

36 

0 

3 

96 

0 

1 

98 

0 

62 

36 

0 

3 

95 

0 

0 

97 

183±7 

2 

82 

16 

8 

7 

83 

9 

4 

86 

1 

82 

15 

2 

7 

88 

2 

3 

92 

325±7.0 

2 

88 

8 

2 

14 

81 

2 

5 

90 

2 

87 

8 

0 

17 

78 

0 

5 

89 

325±9.5 

3 

88 

6 

3 

16 

77 

0 

6 

86 

3 

87 

6 

1 

18 

76 

1 

6 

88 

448±1.5 

0 

17 

81 

0 

3 

94 

0 

0 

96 

0 

14 

83 

0 

3 

93 

0 

0 

96 

448±7.2 

0 

56 

40 

0 

16 

79 

0 

3 

91 

0 

57 

37 

0 

18 

73 

0 

3 

88 

642 

0 

58 

37 

0 

24 

68 

0 

6 

83 

0 

58 

33 

0 

27 

60 

0 

7 

78 

874 

0 

60 

30 

0 

29 

56 

0 

11 

70 

0 

59 

25 

0 

32 

45 

0 

13 

60 
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Figure 1 : The IWP in kg m for the (a) MM5 March 2001 blizzard, (b) the 20 January 2007 lake 
effect snow storm, and (c) the 22 January 2007 synoptic snow event. 
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Figure 2: Vertical profiles of hydro meteor liquid and ice water contents (upper plots) and the 
temperature and relative humidity (lower plots) for the anvil, heavy rain, light rain, heavy snow 

and light snow cases. 
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Figure 3: For this work, the (a) dendrite particle shape used for snow particles, and (b) the six- 
bullet rosette shape used for graupel particles. The non-precipitating cloud ice is assumed 

spherical. (From Liu, 2004) 
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Figure 4: Coefficients for absorption (solid lines) and scattering (dashed lines) for liquid rain, 
dentritic snow, six-bullet rosettes, solid ice spheres, and fluffy ice spheres (10% ice, 90% air) for 
hydrometeor contents of (a) 0.01 g m' 3 , (b) 0.1 g m' 3 , (c) 0.5 g m' 3 , and (d) 0.75 g m°. 
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740 Figure 5: Temperature weighting vectors Wall(z ( )^( z () for the nadir-viewed anvil, heavy rain, 

741 light rain, heavy snow, and light snow cases. Upper panels are for TRMM frequencies (10-89 

742 GHz), middle panels are for AMSU-B frequencies (89-183±7 GHz), and lower panels for the 

743 submillimeter-wave channels (325-874 GHz). 
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Figure 6: Brightness temperature values using a standard clear- air winter atmosphere profile for 
varying surface emissivities and surface temperatures at (a) 10 GFIz, (b) 37 GFIz, (c) 89 GFIz, (d) 
166 GHz, (e) 183±3 GHz, and (f) 183±7 GHz. The X-axis indicates the values used for surface 
emissivity, while the colored lines indicate surface temperature. 
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Figure 7: Percentage contributions from the surface (green), the hydrometeors (red), the relative 
humidity (light blue), the cosmic background (<1% for all channels except 10 GHz where it is 
1.3%)+cloud water+cloud ice (purple), and the atmospheric oxygen and nitrogen (gray) to the 
resultant nadir brightness temperatures for (a) the heavy rain over ocean case, (b) the light rain 

over land case, and (c) the anvil profile case. 
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770 Figure 8: Similar to Fig. 7, percentage contributions for the single profile (a) heavy snow, (b) 

771 moderate snow, and (c) light snow cases. 
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Figure 9: Similar to Fig. 7, percentage contributions averaged over IWP of 0.5 to 1 .0 kg m' 2 for 
the (a) WRF 20 January 2007 lake effect case and (b) the WRF for 22 January 2007. 
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789 Figure 10: Similar to Fig. 7, percentage contributions of a March 2001 snow blizzard MM5 

790 simulation averaged over (a) IWP of 0.5 to 1 .0 kg m' 2 and (b) IWP of 9.5 to 10.0 kg m' 2 . 
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